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Abstract 

A system of partial differential equations describing bilayers am- 
phiphilic membranes is studied by the Lie group analysis. This algo- 
rithmic approach allows us to show all the symmetries of the system, 
to determine all possible symmetry reductions, to recover the axisym- 
metric solutions, obtained by other euristic methods, and finally to 
address the question of new similarity solutions. 

1 Introduction 

In the last decades a great interest has been attracted by the membraneous 
systems in biophysics In particular the binary systems amphiphile- 

water has been considered [^. In water the amphiphilic molecules form sheets 
of costant area, in the form of bilayers, because this minimizes the contact 
of the hydrocarbon tails with water. For a flat piece of membrane, the tails 
are still exposed to water contact at the membranes edge. This leads to an 
effective line tension which, for large enough membranes, pulls the edges to- 
gether and closes the membrane, i.e. we get a vesicle. As reviewed in Sec. 2, 
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a detailed study of such kind of membranes can be performed in a geomet- 
rical setting, following the Helfrich approach. However most of the known 
analytical solutions [|| can be obtained only under particular geometrical re- 
strictions. So the aim of this article is to exploit a systematic group analysis 
(Sec. 3), in order to find the symmetry properties of the system involving the 
shape equation and the Gauss-Codazzi-Mainardi equations. Then, in Sec. 4 
we use the above analysis in performing all the possible symmetry reductions 
and to determine classes of particular solutions. Finally in Sec. 5 we address 
some comments and final remarks. 

2 The Helfrich functional and the vesicles shape 
equations 

The membrane's shape is described by a three-component vector field x(m^, u^) = 
{x^ {u^ , u^) , x"^ {u^ , u^) , {u^ , u^)) depending on the two internal coordinates 
{u^,u^) of the membrane. To be self-consistent, we introduce a bit of the 
geometrical concepts used below The tangent plane to a point P 

on the surface is spanned by two tangent vectors 

ta = SqX, a = 1,2 

where 8^ = 

The two tangent vectors determine the metric tensor or the first funda- 
mental form Vti 

(summation over repeated indices is understood), by which we can express 
the infinitesimal element of distance ds on the surface 

ds^ = Qapdu'^du^ . 

The infinitesimal element of area dS is given by 

dS = a/ det gai3dv}dv? . 
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The metric tensor gap characterizes the intrinsic geometry of the sur- 
face. But for a complete description of its embedding into euclidean three- 
dimensional space we need the second fundamental form 

n2 = bafsdu'^du'^ a,/? =1,2, 

where 
and 

ti X t2 

n 



||ti X tall 

is the unit normal vector on the surface. Embedding of the surface into 
is dictated by the Gauss- Weingarten equations (GWE) 

A ^T<r («,/?,7 = 1,2), (1) 

where g^"' = {g^^)i3-y and F^^ are the so-called Christoffel symbols 
pa <t5 ( ^9as_ dg^ _ dgap_ 

The two most important scalar quantities on the surface are the gaussian 
curvature K and the mean curvature if, given by the determinant and the 
half-trace of the matrix 6^ = g^'^ba-y '■ 

K = det(6^) = -^, (2) 

where Ri and R2 are the principal radii of curvature and they define the 
extrinsic properties of the surface. 

Inspired by the physics of liquid crystals, Helfrich ^ considered the 
derivatives of the unit vector n = (n^,n^,n^) with respect to v} and 
as the curvature strains. The length and energy scales are not the size of the 
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molecular head group, or the molecular energy involved in assembling the 
amphiphilic sheet. Indeed the basic energy is that one needed to bend the 
sheet, i.e. an energy which arises microscopically from a change in the local 
concentration of head and tail groups of the amphiphiles. 

Over short distances, this energy will ensure that the membrane is rela- 
tively flat. However, over larger distances, the effect of thermal fluctuations 
is to change appreciably the orientation of the membrane. The scale of dis- 
tances at which an orientation correlation occurs is called the persistence 
lenght. It can be calculated from the free energy describing a system of 
sheets. 

We assume that the free energy which describes the deformations of the 
membrane depends on its shape only. The local free energy per unit area is 
written as an expansion in powers of local curvatures, and the integrated free 
energy depends parametrically on the unknown coefficients in this expansion. 
These coefficients are the spontaneus, or natural, curvature of the membrane, 
and its bending and saddle-splay moduli. Moreover, the free energy must 
be invariant under translations and rotations and under reparametrization 
transformations of the local coordinates, because of the fluidity assumption. 
Thus, only the invariant linear and quadratic forms can enter into, namely 

dn^ dn^ / dn^ dn^ \ ^ dn^ drP dn^ din? 
du^ dv? ' \ du^ dv? J ' du^ du^ du^ dv} 

Then, the elastic free energy takes the form 

where the so-called bending rigidity k and the gaussian rigidity k are the 
elastic moduli. Making use of the GWE (|l]) and of the expressions (|^-|^) for 
K and H, Fc can be written as 

F^ = k I {H - HsfdS + k[ KdS 
Js Js 

where if^ = — y is called spontaneus curvature and takes account of the 
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asimmetry effect of the membrane in its enviroment. 

When the bending rigidity k ^ ksT, the persistence length is much larger 
than the diameter of a vesicle. Thus we can neglect thermal fluctuations 
and the shape of a vesicle is determined by minimizing the elastic energy. 
Moreover, the area S of the vesicle is fixed for two reasons: a) there is 
very little exchange of molecules between the bilayer and the ambient on 
experimental time scales; b) the elastic energy associated with displacement 
within the membrane surface is negligible on the scale of the bending elastic 
energy. Furthermore, we can assume that the variations of volume, that is of 
the water enclosed in, are negligible. On the contrary, a net transfer of water 
through the membrane would lead to a variation of the osmotic pressure, 
associated with an exchange of the energy in the system, which is huge with 
respect to the bending energy scales. With these costraints, the equilibrium 
shape of a vesicle is determined by the minimization of the Helfrich's shape 
energy F ^ given by 



where A and Ap = pe — Pi are treated as Lagrange multipliers or they can 
be prescribed experimentally by volume or area reservoirs. The membrane 
shapes of equilibrium are defined by the corresponding Euler-Lagrange equa- 
tion (the shape equation) 



where 5*^^-* indicates the first variation of the functional under normal in- 
finitesimal distortion of the surface. The resulting equation is 



Ap - 2XH - 2k{H - HsfH + 2k{H - Hs){2H'^ - K) + kAH = 0, (5) 



where A applied to H represents the Laplace-Beltrami operator on the sur- 



The integrability conditions of the GWE are known from the differential 





= 0, 



face 
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geometry and they are 



- Sf = buTl, + b,,{r\, - r\,) - b,,rl„ (6) 

9^12 db22 _ 7, pi I ^ (-p2 pl ^ h 

'd^ ~ IhIT' ~ "11^ 22 + 22^ ^ 12) ~ "22i- 12- 

These relations are necessary and sufficient conditions for the coefficients Qaf^ 
e bap to be the first and the second fundamental forms of a surface immersed 
in the euclidean space. For this reason they must be included in the set of 
equations, which describe the shape of a vesicle. So the equations (|) and 
(H) describe the vesicle configurations in M^. 

In order to simplify the structure of the system, with no loss of generality, 
we choose metric on the surface in the conformal form 

where (x, y) are conformal local coordinates on the surface. Furthermore, in 
order to put in evidence the role of the mean curvature, it is more usefull to 
introduce the following variables 

1 _ 1 

q = \{bn + b22)e--^e-. (8) 

Thus, the conformal metric takes the form 

ds^ = 4q^ip\dx^ + dy^). 

For sake of clarity, we rename the coefficients bu and bu of the second form 
as 

bu = ^, bi2 = uj, 
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SO the associated matrix is 

' ^ u 
00 8gV(l + Hsip) - ^ 

In conclusion, the system (H-H) takes the form 

q^i^xx + ^yy) + '^q^iQxx + Qyy) " 2v5(g^ + + ^''(Sv? + 02'/'^ + Os'/'^ + a4¥'^) = 0, 

- = {% + ^Lp)q{ipqy + qify), 
uJy + 'd^ = ^q(p{ipqx + q^x) + ^q^qx, 

M^\qxx + qyy) + ^^q^iVxx + ^yy) " 4(^2(^2 ^ ^2) _ 4^2^^2 ^ ^2) _ 
= + _ (8 + £|<^)g2^^9, 

(9) 

where = = 8{2H^ _ |), = 4^ - Si/.f. 

This is a system of nonUnear partial differential equations in two di- 
mensions and in four dependent variables. Its integrability properties are 
unknown: only some particular solutions have been found on the basis of 
euristic considerations (for example by imposing axisymmetry) |[20|,[^,fl. 



3 Group analysis 
3.1 Symmetry algebra 

A specific purpose of the present article is to study the Lie-point symmetry 
group of the system (||). That is, we are looking for a local Lie group G of 
point trasformations on the space of the independent and dependent variables 

X xU = {{x, y, uj, q, (/?) : t.c. : g ^ 0, : ^ 0} C x 

into itself, such that each trasformation g G G acts as 

{x,y,d,uj,q,lp) = g ■ {x,y,T^,u,q,ip) = (Aj, A^, fij, fij), 

where A^ = Ag{x,y,'d,uj,q,Lp), fi^ = Q^g{x,y,'d,uj,q,Lp) are smooth real 
functions, depending also on a suitable set of r real parameters g, conti- 
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nously varying in a neighbourhood of the origin in W. In this section we 



will follow the standard group analysis scheme developped in ||TT|,||T^],||T5|. 
For a symmetry group of a system of differential equations we mean that 
each element g trasforms solutions of the system into other solutions: 
0{x, y),uj{x,y), q{x, y) , ^{x, y) ) is a solution whenever ( y) ,u{x, y) , q{x, y) , v?(x, y)) 
is one. Equivalently, we can say that the symmetry group leaves the system 
invariant. Moreover we are interested in connected local Lie groups of sym- 
metries, leaving aside the problem of the discrete symmetries. Instead of 
considering finite trasformations, in the study of connected groups, we can 
restrict our attention to the associated algebra, that is to the infinitesimal 
trasformations. The symmetry algebra ^ of the system is realized by real 
vector fields, of the form 

ox oy ov ouj oq 0(p 

where ^* and 0q are functions of x,y,d,u,q,Lp. These functions have to be 
determined by the invariance condition of the differential system (P) , namely 

pr^^V [Aj,] = whenever = 0, /or z/ = 1, 4, (10) 

where = indicates each of the equations in (H), and 

a=l a=l y a=l -^-^ a=l -^y a=l Uy 

is the prolonged vectorfield on the space containing the derivatives of the 
dependent variables. Here we have indicated with v} = = uj, = 

q, u"^ = tp and the components 0^ and 0^* can be computed algorithmically 
|[TT|| . The equations ( [T0| ) lead to a large number of linear partial differential 



equations, called determining equations, for the functions C,\ (pa- Solving 
them, we can distinguish two cases: 

I) if the phenomenological parameters are not vanishing {a2, 0.3,04) 7^ 
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(0,0,0), the symmetry algebra Li is spanned by the vector fields 



+C.C., 



where z = x + iy and ^ = .^^ + z^^, being arbitrary real harmonic 

functions satisfying the Cauchy-Riemann conditions 

II) restricting the phenomenological parameters to be (0:2, as, 04) = 
(0,0,0), the symmetry algebra Ln is spanned by the vector fields 



vi(e^(x,y),e'(x,y)) 



V2 



vi(e(^)) = e5.+ 

[-(1? - za;)a^ - (a; + - 4zgV)5c. - \d, 

+C.C., 



+ 



with the same notation used above. The non trivial commutation relations 
for Lj are 

v(e(^)),v(e(^))l =v(ee;-ee.), (n) 



and for Ljj 



vi(e(^)),vi(e(;^))j =vi(ee.-ee.), 

[Vl(e(z)),V2]=0, 



(12) 
(13) 



where ^{z) and ^(2;) are arbitrary complex analitic functions. In both cases 
the symmetry Lie algebras are infinite -dimensional. In particular, the alge- 
bra Li is a realization of the conformal transformations algebra in the plane 



Lc ||18|. Moreover, Ln is the direct sum of Lc and of the one-dimensional 



+ 
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subalgebra Lg generated by V2 



-III 



Lr ® Ln 



A more explicit base for the algebra Lj (Lc) is obtainable by developing the 
functions ^{z) in Laurent series, so that the generic vector field is 



L 6 2. 

+C.C., 



+ 



where a G [0, 27r[. Thus the commutation relations become 

[L„(a), Lm{(i)] = (m - n)Ln+m{a + 



(14) 



In passing, we note that for a = the generators L„(0) span a centerless Vira- 
soro algebra ||ll. The generators {L_i(0), L_i(f ), Lo(0), Lo(f ), Li(0), Lid)} 
span the s/(2, C) algebra, which is the algebra of Mobius transformations in 
the complex plane of local coordinates of the surface. 

In order to give an interpretation of this symmetry structure of the system 
d^), we consider the effect of finite transformations on the metric. The finite 
transformations generated by vector fields v and vion the variables q, z, 
1 are given by 



z{e, z) = T-\e + T{z)), z{e,z) = T-\e + T{z)), 



(15) 



where e is the group parameter and 

r(^) = 

While q and tp transform as 

q{X,q,z,z) = q{z,z) 



ds 



(16) 
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Moreover, being 
we have 



ds^ = Acf'Lp^dzdz = Aq^{z{z), z(z))ip'^{z{'z), z(z))dz(tz = ds^ , 

that is the metric is invariant under group transformations. Moreover, since 
the mean curvature H is related to tp via the eq. (0), by virtue of the eq. ([16D, 
we find that H remains invariant or, equivalently, it transforms hke a scalar 
under the group transformations. Concerning the gaussian curvature it 
is a metric invariant, thus it remains invariant under the symmetry group 
because the metric tensor does. In conclusion we have seen that both H and 
K are invariant under the class obtained of transformations. Consequently, a 
surface, which is a solution of system (H), cannot be "deformed" into another 
one by using a group transformation. 



3.2 Classification of the low dimensional subalgebras 

In order to perform all possible symmetry reductions of @, we must proceed 
to the subalgebras classification under the inner automorphism group (or ad- 
joint group). This leads to build up the so-called optimal system, that is the 
set of the representatives of the s-dimensional subalgebra classes, which are 
pairwise nonconjugated. We start by considering the conjugation classes of 
1-dimensional subalgebras. A first simple result can be achieved by comput- 
ing the adjoint action on a generic element Lm{a) (with m 7^ 0) from the 
commutation relations ( [T^ 



Ad[exp{eLo{a))] (L„(/3)) = e-™^'°'"L^(/3 - m5 sin a), 

with e G M, and a, /3 G [0, 27r[. Then, by a suitable choice of the parameters, 
we can conjugate any subalgebra spanned by {-^^m(/3)}„-^o {-^m(0)}„_^o. 
More generally, this treatment can be applied looking for the action of a 
generic element v(^) on a vector field v(5) (being C = + ^ind 6 = 
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6^ + analytic functions oi z = x \ iy) [0. Integrating the corresponding 
differential equation we are lead to the functional equation 

5(i)=M^(i))- 



where z{z) can be obtained from ([T5|) and 

5(F) |e=o = /u(^) 

for ii{z) given. Now, if we would like to conjugate v(yu(2;)) to the translation 
generator v(l), we have to find a transformation generated by a vectorfield 
v(^), by solving the functional equation for ^: 

1 = M^(2)) 



Now, the rectification theorem for vectorfields ||12] suggests and Neumann 
proved that the real version of this functional equation admits a solution 
for smooth /i and out of a neighbourhood of its zeros. Thus, we can conclude 
that all 1-dimensional subalgebras of L/ are conjugated among themselves, 
at least in a suitable domain. Therefore, we can choose the translation vec- 
torfield L-i(O) = v(l) as the representative of all 1-dimensional subalgebras. 
We now consider a 2-dimensional subalgebra {v(C), v(?7)}. In view of the 
above result, we can conjugate v(C) to v(l) and ^[j]) will be transformed 
into some other v(r^. In addition the commutation relations (11-13) and the 
closure condition require 

[v(l), v(^] = v(^^) = av(l) + fev(rjO, 

that is 



Performing the transformation F= e"''^, we end up with 

= -bzd^, e^'d^ = -bd^. 
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Thus, we can conclude that all the 2-dimensional subalgebras are conjugated 
among themselves and, precisely, to the subalgebra 

Lirf = {L_i(0),Lo(0)}, 

which can be taken as a representative of all 2-dimensional subalgebras. 

3.3 Reduced systems 

In order to find explicit group invariant solutions it is enough to restrict 
ourselves to 1- and 2-dimensional subalgebras. In fact, in these cases we 
have 4 group invariants, which are in one- to- one corrispondence with the 
dependent variables, so that their role can be interchanged. Concerning 
the independent variables, we have respectively 1 and new indipendent 
similarity variables. Thus, the original partial differential system is reduced 
to an ordinary differential one or to a functional one, respectively. 

We now perform the symmetry reductions of (^) starting from the 2- 
dimensional subalgebras, by considering the representative subalgebra Ltd- 
The correspondent subgroup has the following basic invariants 

h = sq, h = tp, 
where s = z — From which we obtain 



= 




h 

Q = — 

s 


_Q 
y 




e 


P 


n 




y 


y2 



where 9, Vt, Q, ipo are arbitrary real costants. By substituting these relations 
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in d^) we obtain 

Q{^o) = ^^oQ' (8 + fv^o), 

n = o, 

The inconsistency of the last equation induces us to conclude that there are 
not group invariant solutions with respect to any 2-dimensional subgroup. 

Now we consider solutions invariant under 1-dimensional subgroups, all 
represented by the translation group generated by the algebra {L_i(0)}. 
With respect to the previous case, we have to add only one more invariant, 
namely s, that is y. Then, the invariant solutions have the form 

1} = i}{y), u = uj{y), q = q{y), = ip{y). 

The reduced system becomes a system of ordinary differential equations 

q'^ + 2q^^ - 2^ (|)' + q' (8y. + a^cp' + a,p' + a,cp') = 0, 

= af + _ (8 + fip) q^p>d, 

where we put to = being a constant. In the next section we will look for 
some particular solutions of this system. 

The importance of this system lies in the fact that all possible 1-dimensional 
symmetry reductions lead to it, modulo transformations of the dependent 
and independent variables. In particular, this has to be true for all systems 
obtained by using some clever parametrization of the axysymmetric surfaces 
studied in the litterature ||g , ||21| , [|| . 

4 Solutions 

In order to find particular solutions of ([T7|), we make some simplifying hy- 
potheses. The aim is to give explicit solutions to the reduced system for the 
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amphiphilic membranes configurations. 



4.1 Sphere 

We first consider the case ip = ipQ. This condition means that we are taking 
under consideration only the constant mean curvature surfaces, since 

^0= rr ^ rr ^ H.^H,. (18) 

-Ho — ns 



The system ([17|) takes the form 



2qip^^ - 2<^o (l)' + (8(/^o + cx^vl + a^^l + a.ipl) = 0, 

fy = + f^o) q^o%, (19) 



The second equation gives 

where B{lpq) = (8 + ^^q) fo- The remaining equations give 

where A{ipo) = {8ipo + a2fl + ctsV'o + <^4V^o)5 with the conditions 

as = 0, 
^0 = 0, 

Assuming A{ipo) ^ 0, the last relation is a condition on the mean curvature 
in terms of the phenomenological parameters, namely 

^R.klil - 2HlkHQ - 2\Hq + = 0. (21) 
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A non trivial solution (modulo translations) of (pO|) is 



'^^^^ C cosh{^/Soy) ' 

where C = 2 \Ho{po\ and 6q is a constant of integration. The gaussian curva- 
ture is 

Thus, the described surface is a sphere of radius R = under the condition 
pl|) . By performing a coordinate change x/ = ^/6oX, yl = ^/Soy, the metric 
and the second fundamental form result 

and the representation in the three-dimensional space is given by 

1 Rsinxf 2/ ^ cos xf n T^ , 

X ixf,yf) = , X [xf,yn = , x ixf,yn = —Rtanhyf. 

cosh yf cosh yf 

4.2 Delaunay's surfaces 

We now consider costant mean curvature surfaces in the particular case of 
Hq = Hg. In this case we cannot use the relation (|18|), so that we are forced 
to perform some slight changes. We first introduce the following variables 



in terms of which the reduced system takes the form 



- (^0 - (f )') + 2^(1 - m (p0 - (I)') + 



4^^P0 + 4pVf - (f ) - 4^ (I) ' = + - 8pV^, 



[22] 
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where (3i = Hs, /32 = f , = -ip- By imposing the condition ip = the 
first equation in (|22D gives 

Ap = 2XHs. 

The second equation becomes 



and the last one takes the form 



P 



(fp ( dp 



dy"^ \dy 



(23) 



where = jf^ii^o + aD- Equation ([2^ ) is a shght generahzation of equation 
([20|) . The integration of (|^) gives 



P{y) = a/t^V 1 - (T'^sn'^{2^/r^{y - yo),cr) 



(24) 



where ro = ' ^ "° , a = / ^VAT^-D ^ r 2|D|a constant of 

integration. Notice that for = 0, that is = cks = 0, we obtain the 
sphere of radius R = jj-. Furthermore, choosing D = iy, we get a cyhnder 
with H = Hg. In the general case, the two fundamental forms are 



ds' 



no 



^[1- ahn' {2^iy - yo),a)] {dx' + dy' 



^[l- ahn'' {2y/¥^{y - yo),(T)] + ^ dx^ + 2a^dxdy+ 



+ 



4r2 
/3 



i[l-a2sn2(27^(y-yo),a)] -^0 dy\ 



Without loss of generality, we can always choose a coordinate system in which 
CKs = 0. Indeed both H and K depend only on 6q and D. In order to simplify 
the last expressions, we perform a coordinate change {x^y) {<P,p) by the 
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1 



Figure 1: Unduloid /9i = 0.5, i^o = 7.2 



equation ([2^) for y and via the relation x = ^0. Thus we obtain 

1 r]^2 I ^2ja2 



'2 ^ .4, 



^2 = ji^-^dp-+ (25) 
These two forms represent a family of surfaces parametrized by -i^o, ^o- In 



the particular case 5o = ^ ^^^^ the two forms (|25D generate the Delaunay 



surfaces p^]||22|. Precisely we have the unduloids when < (31^0 < 1 (fig- 
ure|l]), and the nodoids when /^f-i^o < (figures ||,0,|D. Otherwise, when the 
above constraint on 5q is not satisfied, we guess that a larger class of surfaces 
appears. However we leave this study to the next future. 

4.3 Toroidal surfaces 

Now, we leave the class of the constant mean cuevature surfaces and consider 



the reduction of the system (ITtD under the particular restriction q = Qq 
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Figure 4: Nodoid /3i = 1 = -2.8 
constant. The reduced system (^) takes the form 



f = (8 + f V^) 1^ (26) 



M 
dy 

Avql'i^ - Aql (g)' = al + - (8 + f v^) goV^- 
Integrating the second equation we obtain 

%) = goV (8 + y</^)+ci, 
which substituted in the third equation of pB) gives 



By comparing this equation with the first one of the system (|26|), we obtain 
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Furthermore, the integration of the first equation ( p^) gives 

d(p\ 2 2 ^2 S '-^2 2 4 ^2 5 

So the comparison of the last two equations provides the constraints 

a^ = 0^ Ap = 2XH,, 
f = S ^ A = kHl 
ci = 0, 

= -4- 

Consequently the original system is reduced to the equation 

= -70 - 12^' - 73^^ - = P4(^), (27) 

dy J 872 

2 

where 7o = 72 = Sggi 73 = ^QHsq^. In order to ensure the positivity of 
the r.h.s. of (B^) we must impose the condition 



a^<^(ll + 5v^). 

s 



As suggested by the equation (^7|), performing the coordinate transformation 
{x, y) (x, the two fundamental forms are 



2a., 



^2 = (72 + -j-^^ dx^ ~^ , dxdip + 



70 - 12^'^ - 13^^ - s^V'^ 



4 ( -70 - 72V?^ - 73^^^ - s^V^^ 



with 0</?<(/J<aifi?s<0 and (3 < Lp < a < if Hs > 0, being a, /3 the 
two real roots of ^(v?), which is positive in jo;, The difference in the sign 
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of Hs depends only on the different orientation of the normal vector n. The 
gaussian and mean curvatures are given by 



V 272 



K 



270 , 27; 



+ 



+ 



71 



72 72 V' 47; 



2 • 



In the special case = and rescaling the coordinates x = -^0, ip = V2p, 
we have 



H 
K 



dp" 



, (28) 

( H,p' + V2p) de' , (29) 

(30) 

m + ^^. (31) 



2p ' 
r2 . V2Hs 



If we interpret the symmetry variable 9 as the rotation angle through the x^- 
axis and p as the radius a/ (x^)^ + (x^)^, then ( p8| - p9| ) are the two fundamental 
forms of the so called Clifford torus, with cartesian equation 



{x'Y + (xY + {x'Y + 



[{x'f + ix^] . 



This vesicle configuration has been detected experimentally by M. Mutz and 
Bensimon |p4| , ||2S| (see figure ||). 



4.4 Circular Biconcave Discoid 

In the previous sections we have considered the special cases in which the 
variables q 01 (p are taken to be constant. In this section we will consider 
the general case in which both q and ip are not constant functions, but we 
suppose that q is functionally depending on ip. Thus we put without loss of 
generality 

. = .M = 4^. (32) 
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Figure 5: Clifford torus 



where p is a function to be determined. We now look at the system ([17|) 
and consider the second equation. Assuming a one-to-one correspondence 
between p and y?, we have 

where the integral is still an undetermined function of p. However one could 
play with in order to get closed expressions. For example, the toroidal case 
can be recovered by assuming tp and p to be proportional, or the expression 
under integral of the form ^. Another interesting possibility arises making 
the following hypothesis: 



1 dip c 
ip'^ dp p' 



(33) 



or, equivalently 



Consequently we have 



cln 



PO J 
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If we now substitute the two last expressions in the remaining equations of 
([T7|), we obtain the compatibility conditions 



a2 = 12c, as = 4c^, 04 = 0, = 'do = 0, 



(34) 



thus c = 2Hs and the other phenomenological parameters A and Ap are 
equal to zero because of (0). Furthermore, we have the following differential 
equation for p{y) 



d 



dp " 
dy 



- +cVln - 

dy \ p \po 



p 



l + ln( ^ 

Po 



0. 



(35) 



In order to integrate this equation we perform the hodograph transformation 

dy 



dp 



rHp)- 



(36) 



The equation (|35|) provides a differential equation for f{p) 

P 



4f _ 2 
dp p 



_ - / = _2cVln' ^ 



1 + ln 
PoJ V VPo 



the general solution of which is 



2„2i 2 / P 



/(p) = pM/o-cVln^(^ 

Po 



{foe 



Thus, the problem of the system (ITTD with tha assumptions (|32| - p3D is re 



duced to the quadrature of (|36D . Furthermore, we can conclude that the 
solution surface is defined by the following two fundamental forms in the 
local coordinate (x, p) 



ds^ = p^dx^ + -^-—dp^, 



fipY 



p 



p 



dp\ 
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Figure 6: Circular Biconcave Surface 



From the above expression of the metric, we interpret p as a radius and x as an 
angular variable, obtaining an axisymmetric surface. A general parametriza- 
tion of this kind of surfaces is given in terms of angle iJj between the tangent 
to the contour and the radial direction 0| 

x^ = cos(x) = sin(a;) 

x^ = J ta.n'il){p)dp + zq 

where sec^ i^ip) = The required integration for x^, which can be numer- 
ically performed, leads to the well known |^ biconcave discoid surface (see 
figure ^ where we set /o = 1 in order to avoid a conical point in p = 0. This 
surface may represent a static configuration for the red blood cells. 

5 Conclusions 

In this article we have shown that the symmetry algebra of the system (|) is 
the algebra of the conformal transformations in two dimensions, that is they 
generate arbitrary reparametrizations of the local coordinates {x,y), accom- 
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pained by the appropriate transformation of the dependent variables. In the 
particular case Hg = Ap = A = the functional F (|^) becomes the so-called 
Willmore functional and the symmetries include a one-dimensional algebra 
Lg (spanned by the vectorfield V2), which generates the scaling transforma- 
tions in the ambient space One of the main result of our research is 
that in the general case (the most important from the physical point of view) 
the group analysis leads to a unique symmetry reduction, that is the system 
([TTI) . By making different simplifying hypothesis (g = qo, ip = ipo, etc..) we 
have recovered the known solutions of toroidal and constant mean curvature 
type, respectively. However, the first and second fundamental forms contain 
possible generalizations of these surface solutions. Moreover we have seen 
how to obtain in a systematic way other classes of solutions, by imposing 
particular relations between q and tp (see eq. p^-|33|). In this way we have ob- 
tained the so-called circular biconcave surfaces, but we presume to be able to 
reduce to the quadratures many other classes of surfaces in the next future. 
Besides, in order to find new solutions and study their stability, we plan to 
use the discrete symmetries of the system, numerical methods in determining 
solutions of the GWE (|I]) and finally we plan to study the second variation 
of the functional F on the solutions. 
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